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We compute statistical properties of the lowest-order multipole coefficients of the magnetic field 
generated by a dynamo of arbitrary shape. To this end we expand the field in a complete biorthogonal 
set of base functions, viz. B — a*" {t)b'' (r) . The properties of these biorthogonal function sets 
are treated in detail. We consider a linear problem and the statistical properties of the fluid flow are 
supposed to be given. The turbulent convection may have an arbitrary distribution of spatial scales. 
The time evolution of the expansion coefficients a* is governed by a stochastic difi^erential equation 
from which we infer their averages {a'^), autocorrelation functions {a^ {t)a^* {t + r)) , and an equation 
for the cross correlations {a^a^*). The eigenfunctions of the dynamo equation (with eigenvalues 
Afe) turn out to be a preferred set in terms of which our results assume their simplest form. The 
magnetic field of the dynamo is shown to consist of transiently excited eigenmodes whose frequency 
and coherence time is given by 3Afc and — l/5RAfc, respectively. The relative r.m.s. excitation level 
of the eigenmodes, and hence the distribution of magnetic energy over spatial scales, is determined 
by linear theory. An expression is derived for P)/(|t[''P) in case the fundamental mode 6° has 
a dominant amplitude, and we outline how this expression may be evaluated. It is estimated that 
(|a'°|'^)/(|a''p) ~ where A'' is the number of convective cells in the dynamo. We show that the 
old problem of a short correlation time (or FOSA) has been partially eliminated. Finally we prove 
that for a simple statistically steady dynamo with finite resistivity all eigenvalues obey KAfc < 0. 
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1. INTRODUCTION 

The origin of the magnetic field of the Earth and the 
Sun is well understood at the qualitative level. Helical 
convection acting on the toroidal component of the field 
generates new poloidal field, while the toroidal field is 
regenerated either by shear flows acting on the poloidal 
field or by the same helical convection now acting on the 
poloidal field. These dynamo processes are in principle 
able to balance resistive decay and to maintain a mag- 
netic field for very long times. The large scale magnetic 
fields observed in galactic discs are likewise believed to 
be due to similar dynamo processes Self-consistent 
hydromagnetic simulations that became available since 
1995 have confirmed this dynamo picture for the geomag- 
netic field 12]. Many groups have since then published 
numerical geodynamo models. Even though the avail- 
able computational means do not permit the parameters 
of the models to be 'earth-like', the magnetic field of the 
simulations has many properties in common with the ob- 
served geomagnetic field 0, H, S H @| • For the Sun with 
its much higher magnetic Reynolds number such simula- 
tions are not yet feasible. 

The availability of numerical geodynamo models has 
opened up the possibility for a detailed diagnostics of 
dynamo action. Kageyama and Sato [7] and Olson et al. 
[ii found that the regeneration of poloidal and toroidal 
field resembles the a^-dynamo scenario of mean field the- 
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ory 0. Wicht and Olson [T^] analysed the sequence of 
events leading to a reversal in a simple numerical dy- 
namo model. Schrinner et al. (llj have inferred mean 
field tensors aik and PiM and the mean field {B) from 
simulations. The mean field is then compared with the 
mean field predicted by the dynamo equation. Reason- 
able agreement was found for a simple magnetoconvec- 
tion model and for simple dynamo models. For a recent 
review we refer to Wicht et al. p^ . 

The evolution of the magnetic field B in the conduct- 
ing fiuid of a dynamo is governed by the induction equa- 
tion: 

_ = Vx(y-77V)xB. (1) 

The dynamo is located in a volume V with exterior vac- 
uum E, see Fig. [TJ The shape of the dynamo need not 
be spherical. The flow V consists of a stationary compo- 
nent V and the turbulent convection u, which may have 
an arbitrary distribution of spatial scales: 

V = v + u . (2) 

A popular line of attack is to average over the turbu- 
lent convection u and to derive an equation for the mean 
field {B) (the so-called mean field dynamo equation or 
briefly dynamo equation [3]). We follow a different 
path, and we expand the field B in a complete set of 
functions b^{r). We then determine the statistical prop- 
erties of the expansion coefficients. The mean field (B) 
will appear only occasionally, as a mathematical concept 
without much physical meaning attached to it. 
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The idea to study dynamo-generated magnetic fields 
by an expansion in multipoles goes back to Elsasser [3] . 
We follow the same technique and obtain a set of equa- 
tions for the expansion or mode coeSicients. These equa- 
tions contain a random element as the fluid motion V 
in the induction equation consists of a steady part with 
a superposed turbulent convective component. The new 
aspect is that we use the theory of stochastic differen- 
tial equations [l5j to infer the statistical properties of 
the mode coefficients. We consider a statistically steady, 
saturated dynamo with a selfconsistent mean flow u and 
turbulent flow v. 

We treat a linear problem and consider u and v as 
given. Until recently it had been tacidly assumed that 
the solution of the induction equation ([1]) represents the 
selfconsistent field B obtained from a nonlinear solution 
of the MHD equation (provided one uses the exact self- 
consistent flow u + v). But we now know that this is 
not correct. There are statistically steady, saturated dy- 
namos whose velocity field, taken as a given input fiow in 
the induction equation, acts as a kinematic dynamo with 
exponentially growing solutions 17?|. In those cases 
the induction equation on its own is obviously unable 
to reproduce the selfconsistent field B. Many questions 
regarding this unexpected phenomenon remain to be an- 
swered. For example, it is not known to what extent it is 
a universal feature. Schrinner and coworkers (to be sub- 
mitted) have found several counterexamples, and their 
results suggest that the flow fields of fast rotating geody- 
namo models are also kinematically stable. The solution 
of the induction equation with the (selfconsistent) fiow 
taken from these dynamos is (after a transitory period) 
up to a constant factor equal to the selfconsistent field B, 
independent of the initial condition. In the absence of a 
generally agreed-upon terminology we shall in this paper 
refer to these dynamos as kinematically stable dynamos. 

Here we restrict ourselves to kinematically stable dy- 
namos, so that the solution of H]) faithfully represents the 
actual field B. Otherwise, the dynamo model is general 
and may be of the geodynamo or solar type. We take the 
existence of a linear dynamo instability and its nonlinear 
saturation for granted, and study the fiuctuations of the 
system around this nonlinear equilibrium state driven by 
the turbulent convection. A proviso must be made for 
galactic dynamos since these may not yet be in a statis- 
tically steady saturated state. It is therefore not clear if 
the theory developed here is applicable to these dynamos. 

In this study we do not focus on the mean field con- 
cept (B); we rather expand the field B in a complete 
set of functions b^{r). But as we determine the statisti- 
cal properties of the expansion coefficients, notions from 
mean field theory pop up simply because we compute av- 
erages. For example, the dynamo parameters a and f3, 
well known from mean field theory, emerge because they 
are connected to the simplest nontrivial averages of the 
turbulent convection. It is almost unavoidable that they 
should appear in any theory that considers averages over 
the turbulent convection. In other words, as we compute 



the statistical properties of the expansion coefficients we 
make contact with mean field theory, but we do not use 
the mean field concept (B), except occasionally in pass- 
ing. 

The context of this study is as follows. Hoyng et al. 
[m analysed the statistical properties of the multipole 
coefficients of the mean field of a simple mean field geo- 
dynamo model excited by fluctuations in the dynamo pa- 
rameter a. Here we addess a much more general problem: 
the statistical properties of the multipole coefficients of 
the field itself due to forcing by the turbulent convection. 
The present paper builds on an earlier study of Hoyng 
[TqI that addressed the same problems, but with lim- 
ited success. The question of the r.m.s. mode excitation 
level was only solved conceptually, and the derivation 
of the time evolution of the mean expansion coefficients 
contained errors. Moreover, the study was restricted to 
dynamos with homogeneous isotropic turbulence. These 
points of criticism have been removed here. Since we 
study the relative excitation level of magnetic overtones 
in the dynamo, the present work may also be regarded 
as a spectral theory [20| yielding the distribution of mag- 
netic energy in a finite dynamo as a function of spatial 
scale (i.e. mode number), with due allowance for the 
boundary conditions. 

This paper is organised as follows. In Sec. Owe summa- 
rize the properties of the function sets we use to represent 
the magnetic field of the dynamo, in particular biorthogo- 
nality and closure relations. Then we derive the equation 
for the mode expansion coefficients in Sec. [31 Next we 
determine the statistical properties of the expansion co- 
efficients. In Sec. [5] we derive the time evolution of the 
mean mode coefficients and show that it is governed by 
the operator of the dynamo equation. In Sec. [S] and [7] we 
compute the excitation level of the overtones relative to 
the fundamental (dipole) mode. Finally, in Sec.[8]we dis- 
cuss our results and potential applications in a somewhat 
wider perspective. 

2. COMPLETE FUNCTION SETS 

We summarize here some technicalities concerning the 
function sets that we emplo y. Some of these properties 
have been derived elsewhere [l9l . [2l| , and we collect them 
here in the interest of a homogeneous notation, taking 
the opportunity to correct some mistakes. Readers may 
move directly to the next Section, and return here for 
reference. 

We expand the field B of the dynamo, see Fig. (TJ in a 
complete set of functions tr' (r) : 

B = Y,a^t)y{r) . (3) 

3 

These functions are often eigenfunctions of some differ- 
ential operator, and since this operator is usually not 
self-adjoint the functions are not orthogonal. This prob- 
lem is handled by using the adjoint set b*(r), with the 
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FIG. 1: A dynamo in a volume V that need not be spheri- 
cal, with exterior vacuum E. The fluid in V has a magnetic 
resistivity rj, and the flow V consists of a stationary compo- 
nent V and the turbulent convection u, which may have an 
arbitrary distribution over spatial scales. The currents J are 
restricted to V, and both currents and flows are tangential 
to the boundary dV. The magnetic field B is continuous 
through dV. 



following property 

/ b' lt> d^r = S'^ . (4) 

Jv+E 

The base functions b*(r) and b-' (r) constitute a 
biorthonormal set [l^, ; they are continuous through 
the boundary dV of V, just Hke B, and potential in E. 
A few remarks on the notation: 

- the hat " includes complex conjugation, so 

13 = ^a"(i)b^(r), 

i 

- upper indices enumerate the base functions, 

- lower indices indicate vector components, 

- currents are defined as J = V x S - the factor An/c is 
absorbed in J. 

Examples of such function sets are the free decay 
modes, or the eigenfunctions of a homogeneous spheri- 
cal dynamo (see Ref. [oj. Chap. 14). These two sets are 
actually self-adjoint so that the adjoint sets are obtained 
by complex conjugation: 6*(r) = b"(r). Usually the 
situation is more complicated, and we refer to Schrin- 
ner et al. [1^ for the construction of the adjoint set of 
the eigenfunctions of the dynamo equation. From that 
construction a number of properties emerge, such as j = 
adjoint of V x 6 = V x b, i.e. the adjoint operation 
commutes with V. 



2.1. Biorthogonal sets in V 

It is of course possible to define an inner product in 
V alone, and to construct an adjoint set of so that 
these are biorthogonal in V: Jyb^ ■ b-' d^r = 6^^. But 



as we shall see these sets are not very convenient. Fortu- 
nately, the boundary conditions permit construction of 
a very useful biorthogonal set in V, starting from sets 
{b'} and {b-*} that are biorthogonal in V + E. For two 
fields bi and b2 with associated currents j = V x b 
and vector potentials a, where V x a = b, we have 
V • bi X a2 = ji • a2 — bi • b2, whence 

/ bib2d3r = 

Jv+E 

/ ji-a2d^r - (* +* )bixa2-dV. (5) 

Jv \JdE JdvJ 

Here dV and dE indicate the boundary of V and E, 
respectively. The surface integrals cancel because the 
field is continuous through dV and the vector potential 
can be constructed to be so. Volume integrals containing 
currents may be restricted to V since j = in E. This 
leads to: 

/ b* • b^ d^r ^ [ y ■ 

Jv+E Jv 

= [ a'- f dV = 5'^ . (6) 
Jv 

These relations are invariant under a gauge transfor- 
mation a ^ a + Wtp because Jy j ■ Wtp d'^r — /^^ V • 
{4>j)d'^r — JgyipJ ■ d^cr = 0, as (adjoint) currents run 
parallel to the boundary. Currents and vector potentials 
thus form a biorthogonal set in V, and this is essential 
for our study. 

2.2. Closure relations 

Completeness means that any reasonable field B 
should have a unique expansion ([3]). Completeness 
implies the existence of closure relations. The best 
known example is from quantum mechanics where or- 
thonormality of two eigenstates m and n means that 
J tpmiryipnir) d^r = Smn, and the closure relation is 
'4'n{i')*^n{i'') = S{r — r') provided the sum is over 
a complete set of eigenstates 24] . In the present context 
it can be shown that 

^fc2 ^ei — ^ fk2 ^\\ = ^ j\\ 

i i i 

= &{1 -2)5u ■ (7) 

Notation: 1 = r^, 2 = r2, b^^ — £-th vector component 
of b*(ri), etc. Relation ([7]) says that the expressions on 
the left are unit operators in the discrete and continuous 
indices. There are several other closure relations, and the 
one we actually need is 

J23l2b],i = ekimSil ~2)\/r^2 ■ (8) 
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Relations ([7]) and ([8]) hold under an integral sign 
Jyd2--- and are to be regarded as an operator acting 
on all r2 dependence in the integrand. Care is needed 
when applying these closure relations, and for details we 
refer the reader to Appendix C of Ref. [2l| . Recall that 
we have absorbed the factor 47r/c in the definition of the 
current. 



3. MODE EQUATIONS 

On taking the inner product of ([3]) with b*^ we obtain 
with the help of j!]): 



a'^it) = ■Bd'^r 

'v+E 



(9) 



To obtain the evolution equation of a*^, we take the time 
derivative ' = dt and split the integral in two parts: 



B d^i 



(10) 



In the first term we may use the induction equation, but 
we have no equation for dtB in E, where B is simply 
the potential field continuation of B on dV. Omitting 
the second term is no option as it amounts to saying that 
dtB = in E. A way out of this dilemma is to integrate 
the right hand side of ^ by parts, as in ([5]) , to obtain 
(V X A = B): 



a\t) 



f ■ A d^r 



(11) 



Derivation of the mode equations is now straightforward. 
Take the time derivative of and insert the uncurled 
induction equation, dtA — {{v + u) — 77V} x B + Vif}, 
to obtain 



{{v + u) - 77V} X B dV . (12) 



The gradient term produces zero, see below At this 
point we note that if we had used a biorthogonal set 
6'(r), 6^ (r) in the dilemma mentioned a moment ago 
would not occur, because the integrals are restricted to 
V . However (1121) would then read 



a^{t) = / b'^' • V X {(v + m) - 77V} X B dV . (13) 
Jv 

There is no conflict between ([T^ and because they 
use different bases: biorthogonal vnV + E and biorthogo- 
nal in respectively. So 6'^, 2^ in are not the same 
as 6^ f in Uni)! 

The extra V on the right hand side of is very 
inconvenient and cannot be removed by integration by 
parts because the surface integral over dV is nonzero (see 
Appendix |^ . It would stay with us and cause havoc in 
all subsequent analytical and numerical work. Therefore 



((12)) is much to be preferred. This point was one of the 
main drivers for developing the formalism of the previous 
Section. Because of its importance, we discuss the issue 
in more detail in Appendix [X] 

On inserting B — J^e'^'^^^ (fT2)) . we arrive at the 
mode equations 



d'= = J2{S''^ + C''\t)}a'^ , (14) 
e 

S''^ ^ f 3''-{v- T]V)x d^r ; (15) 
Jv 



with 



C'^^(t) 



■ ux d^r 



(16) 



u-j^xb^ d^r 



The base functions that occur in the definition of the ma- 
trix elements S^^ and C^^ may be thought of as spatial 
filters with positive and negative signs indicating how 
the velocities v and u are to be weighted in the vol- 
ume integration over V . The boundary conditions make 
that the integrals are of the type (current|operator|field) 
rather than (fieldj operator] field). The boundary condi- 
tions have therefore a profound effect on the definition 
of these filters. Mode equations of the type (fH|) had 
already been derived by Elsasser jl3|, although he con- 
sidered steady flows only. Compared to the time scales in 
the steady matrix S^^ , the elements C^^ fluctuate rapidly 
because u{r,t) does. The mode equations are therefore 
stochastic equations, and we shall use the theory of these 
equations to infer averages of the mode coefficients. 



STOCHASTIC EQUATIONS WITH 
MULTIPLICATIVE NOISE 



Consider a linear equation of the type 



dx 
'dt 



{A + F{t)}x 



(17) 



where linear means linear in x. The quantity x may be 
anything, a scalar, a vector such as B, or other. The 
operator A is independent of time but F{t) fluctuates 
randomly, with a correlation time Tc and zero mean, 
(F) — 0. This 'noise term' F{t) is said to be mul- 
tiplicative as it multiplies the independent variable x, 
and that makes (fTT]) more complicated than for example 
dx/dt = Ax + F{t), where the noise term is additive. 

A tenet of the theory of stochastic differential equa- 
tions is that the average (x) obeys a closed equation [1^ : 

^(.t) = |a +^"dr(i^(i)exp(AT). 

■F{t-T))cxp{-AT)]{x) , (18) 
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The right hand side of (|18p is actually the first term of a 
series with terms of order (_Ftc)" containing higher-than- 
second-order correlation functions [l^. Here F stands 
for -Fr.m.s. • We assume that the correlation time is small, 
Ftc <C 1, so that only the first term survives and only 
second order correlations matter, as in (fT8|) . Frequently 
also Atc <ti 1 holds, and then we may ignore the ex- 
ponential operators in (|18p . These evolution operators 
exp{zLAT), incidentally, ensure that the transport coef- 
ficients for (a;) derived from are Galilean invariant 
[26t . Relation (fT5|) will be applied several times in what 
follows. 



5. THE AVERAGE MODE COEFFICIENTS 

The equation for the average (a*^) is in principle avail- 
able in the literature, but hidden in two papers [l^, |2l| . 
We summarise the derivation here, and apply (|18p to 



(HI: 



{a^) = I -5^^ + y dr {C'''{t)C'\t - r)) | (a^ . 

(19) 

Summation over double upper indices (here i and £) is 
implied. Before we elaborate the correlation function we 
verify the conditions to be satisfied: 

(1) C^^Tc ^ 1, the usual condition of a short correla- 
tion time, now in the context of the present application 
[C^^ is understood to indicate a characteristic magni- 
tude]. We estimate the value of C'^^t^ for the geody- 
namo. The volume integration in reduces u effec- 
tively to Wr.m.s.^"^^^ where N ^ AR^/\l ~ number of 
convection cells; R is the radius of the outer core, and 
Ac the correlation length (the volume of the inner core is 
neghgible). We estimate ~ V x ^ /{R/k) and 

h'^ b^(^r = 0(1), and Ur.m.s.Tc ~ Ac. This leads to 



-^r.m.s. 'c 



"aT 



-1/2 , yk f^i 



R/k 



d^r 



k fx 



R 



5/2 



(20) 



For the geodynamo we have i?/Ac ~ 3.5. It follows that 
C^^Tc ~ 0.02A: so that C^^Tc <C 1 seems well satis- 
fied for the lowest multipole coefficients a'^ (even though 
Wr.m.sTc ~ Ac !). There is an interesting connection here 
with the problem of the First Order Smoothing Approx- 
imation (FOSA) in mean field theory, and we refer to 
Sec. [S] where we illustrate how this FOSA problem seems 
to be partially eliminated. 

(2) S^^T^ < 1. The magnitude of S*" is set by the 
resistive term in (jlSp . and in much the same way we 
obtain 



For the geodynamo we have 77 ~ 1 m^ s""'^, and S'^^Tc ~ 
3 • lQ~^kl so that S'^^'Tc <€. 1 holds except for high order 
modes. We have to assume therefore that these modes 
contribute little, which is plausible because they will have 
small amplitudes. 

The correlation function in pQ]) is evaluated in Ap- 
pendix [B] and as a result we may write pop as follows: 



kli 



with 



jjki ^ I f . Db' dV , 



(22) 



(23) 



where the operator D is defined as 

Db = vxb + ab-f3: (V6) - 77V x b . (24) 

Notation: {a ■ b)i = a^fefofc and [/3 : (Vb)]i = Pikd'^kbe). 
The operator D is intimately related to the dynamo equa- 
tion, which reads f Appendix [B|) : 



I (S) = V X D{B) 



(25) 



The fact that the operator of the dynamo equation for 
the mean field emerges in (I22p is quite independent of 
which set of base functions is used. At this point we see 
how mean field theory emerges as we study the statistical 
properties of the mean mode coefficients. It happens be- 
cause we take averages, and since the simplest averages 
over the turbulent convection involve v, a and /3, it is 
perhaps not surprising that the operator of the dynamo 
equation for the mean field appears. But we use the mean 
field (B) as a mathematical concept only, without much 
of a physical meaning attached to it. 



5.1. Dynamo action and dynamo equation 

The fact that the dynamo equation for the mean 
emerges as we set up the theory indicates a deep con- 
nection between dynamo action in flows with a random 
element and the mean-field dynamo equation. What this 
relation is can be seen by adopting the eigenfunctions of 
the dynamo equation as the complete function set {b*}, 
with eigenvalues Xi |27i]. Then we have 



V X Db' = A,- b' 



kiriTc/R^ 



(21) 



(26) 



Relations through (|^^ below hold only if we use the 
eigenfunctions of the dynamo equation as our function 
set. Schrinner et al. [2^ show how these eigenfunctions 
and their adjoints may be constructed. 

The computation of D''^ from ([^5]) has several potential 
pitfalls. Integration by parts to Jyb'^-Vx Db^ d'^r, with 
hopes of using (pS]) . is not possible because the surface 
term tf^y b'^ x Db ■ (Per does not vanish. And to argue 
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that ^ = Jy^j^ f ■ Db^ d^r = Jy^j^ b'^ • V x Db'^ 
does not work either because D is not defined in _E. A 
safe way is to uncurl ([26|l to Db^ = Ai a* + V0, and to 
insert that in 

D'^^ ^ [ f ■ d^r = \i 8^^ . (27) 
Jv 

The gradient term produces zero, see below The 
detour via the vector potential in ((27)) allows us to avoid 
trouble due to the discontinuity of the operator D on dV . 

Returning to the evolution of the mean mode coeffi- 
cients, we now have 

= Afc(a'=) , or (a'^') cx exp(AfeO . (28) 

At 

If we choose the eigenfunctions of the dynamo equation 
(the operator V x Z?) as a basis to represent the magnetic 
field, then the mean mode coefficients decay as simple 
exponentials at a rate set by the eigenvalue of the mode. 
Representation of the field B of the dynamo in terms of 
the eigenfunctions of its dynamo equations is therefore 
an efficient way to expose the underlying physics. The 
autocorrelation of the expansion coefficients becomes: 

{a\t + T)a!^*{t)) = .exp(AfcT) (29) 

(no summation over fc). The proof is simple and not given 
here. The mean squares of the expansion coefficients are 
computed in Sec. [71 relation ([iU)) . 

Of course one may always use a different representa- 
tion. In that case D^^ is no longer diagonal and the mean 
mode coefficients have a more complicated time depen- 
dence. Equation may be integrated to 

(a'=)|t = [exp(i?i)]"(«'>|t=o . (30) 

The time dependence of {a}^)\t and the autocorrelation 
function is now a mixture of exponentials. 

5.2. Absence of growing modes 

For the theory developed here to make sense it is nec- 
essary that all eigenvalues have a negative real part, 
3?Afc < 0, and we shall assume that this is the case. 
The implication is that if we measure the flow in a dy- 
namo (model) with all feedbacks of the Lorentz force act- 
ing on it, we should find, on solving (|26p with onj and 
/3ijfe and mean flow computed from these measurements, 
that all eigenvalues have SRAfe < 0. This is one of the 
checks on a correct determination of the dynamo coef- 
ficients. Of course this can only be true if we restrict 
ourselves to kinematically stable dynamos (see Sec. [T] for 
a definition). A proof that 3?A/c < has beengiven for a 
few special cases with infinite conductivity [28l . |29| , but 
a general proof is lacking. In Appendix [Cl we push the is- 
sue one step forward and prove that all eigenvalues of the 
dynamo equation of a statistically steady dynamo with 



locally isotropic turbulence and finite conductivity have 
negative real parts. 

A consequence is that all mean mode coefficients ap- 
proach zero, (a'^) | and that the mean field (E) = 
will also be zero, 

(B) = . (31) 

Ultimately, this may be seen as a consequence of the fact 
that both B and —B are solutions of the MHD equa- 
tions for given V and since there is always a (possibly 
very small) transition probability between the two states, 
a full ensemble average {B) of a statistically steady dy- 
namo will be zero. This is also true if (•) is a time average, 
provided the average is over a sufficiently long time, so 
that all dynamical states of the system are sampled. For 
the geodynamo that would mean an average over a time 
interval containing a large number of reversals (or cycle 
periods in case of the solar dynamo). 

From a practical point of view, it is therefore not ob- 
vious what a time- or an ensemble-averaged field {B) 
is telling us about the actual field B of a dynamo, and 
this is one of the main motivations for the approach pre- 
sented in this paper. A similar objection may be held 
against the mean mode coefficients (a'^). They also do 
not contain much information since they are zero (after 
some time) , but they are indispensable for setting up the 
theory. The cross correlation coefficients {a^a^) that we 
compute in the next Section are much more useful in this 
regard. 



5.3. Mode excitation and phase mixing 

The interpretation of these results is the following [1^ . 
The eigenmodes of the dynamo equation serve as a kind 
of normal modes of the dynamo. These normal modes 
are transiently excited and have a frequency 3Afc. The 
relative r.m.s. mode amplitudes will be computed in the 
next Section. Random phase shifts due to the stochastic 
forcing render the modes quasi-periodic with a coher- 
ence time — l/3fiAfc. Due to the averaging, these phase 
shifts appear as a damping of the mean mode amplitude 
(a'^). This is a familiar phenomenon in statistical physics 
known as phase mixing. 



6. MEAN SQUARE MODE COEFFICIENTS 

In this Section we are interested in the cross-correlation 
coefficients {a'^a^*). The strategy is to derive an equation 
of the type (d/dt) {a'^a^*} = operator - {a'a^*}. We may 
then apply relation to find the equation for {a^a^*) 
(see Ref. ^5|], § 16). The analysis is straightforward (we 
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use ■ = d/dt ): 

+ C^^S^C^^^^) a'a^* . (32) 



We see that the dynamo operator D also plays a role 
in the evolution of the mean square mode coefhcients. 
Relation l|33p becomes more transparent if we adopt the 
eigenfunctions of V x D as our basis: 

(^j^kmin _^ ^tnkm*^ (a^o"*) . (34) 

(summation over n, m, not over /c, 1). The matrix 7Vf*="'^" 
is defined as 

/"OO 

Mkmin ^ dr(C'='"(i)C^"*(i-r)) . (35) 
Jo 

6.1. Interpretation of Eq. ([34]) 

Equation ([34]) describes the time evolution of the mean 
square mode coefficients {a^a^*). We replace d/dt by A, 
as usual, to obtain the eigenvalue problem (the double 
indices fc, £ and m, n may each be grouped into a sin- 
gle new index). Of interest is the largest eigenvalue A, 
which should be approximately zero, otherwise {a^a^*) 
would either grow indefinitely, or become zero. This is 
again a matter of using the correct values of the dynamo 
coefficients. For instance, in a statistically steady geo- 
dynamo simulation the magnetic energy is on average 



This is of the form ((TT]) with x = {a^a^*}. The condi- 
tions of a short correlation time Ft^ <C 1 and At^ <^ 1 
are the same as (f20| and (f2T|) , respectively. The applica- 
tion of recipe for the average is mainly a matter of 
bookkeeping of upper indices: 



(33) 

I 

constant, which implies that the mean square mode co- 
efficients (a'^a^*) are constant, even though some of the 
modes may be quasi-periodic. A determination of the 
dynamo coefficients and Pijk from a measurement of 
velocity correlation functions in this dynamo must then 
result in Eq. having a stationary solution. This solu- 
tion is the eigenvector belonging to the eigenvalue A = 
and it specifies the relative magnitudes of the correlation 
coefficients {a^a^*) (the absolute level is out of reach and 
requires explicit inclusion of nonlinear effects). These 
contain a lot of information, such as the distribution of 
magnetic energy over spatial scales. All other eigenval- 
ues A should have negative real parts, and correspond to 
transient initial states. 



6.2. Computation of m'""'^" 

Expression ([55]) for M'^™'^" looks very similar to the 
correlation function in (|19p . but there is an important 
difference: there is no internal summation over modes. 
The closure relation is therefore no longer a comrade in 
arms, which renders evaluation of psp much less straight- 



dt 



■ (^C'"\t ~ t) ,5^" + C^"* {t - r) <5*™) ) I (a™a"*) 

.kraxin _^ gin* gkm _^ j . _ ^^^km ^ ^^^^^ ^ _ ^^yn* ^krn^j 

dr ((C'='"(t)C^"*(i-T)) -I- (C^"*(t)C'^''"(t-r))U(a™a"*) 



jjkmgen ^ jjen*gkm + dr ((C'=™(t)C^"*(t - t)) + (C^"*(t)C*^'"(i - r))) | (a^a"*) 
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FIG. 2: Intuitive picture of excitation of overtones in case of 
a dominant fundamental mode. 



forward. We insert (1161) in (l35l): 



M 



km£n 



/ dT(C'=™(OC^"*(t-T)) 

Jo 



92 



dT«i<r)- (36) 



Here 1 — ri, 2 = r2, index pi indicates p-th vector 
component at position ri, u^i — p-th vector component 
of u at position ri and time t, etc. Summation over the 
vector indices p, q is implied. 

There is an opportunity to simphfy relation (|36p if the 
correlation length Ac of the flow u is much smaller than 
the size of the dynamo - in other words, if there is a 
separation of spatial scales. In that case we put r2 — 
Vi+p and note that (UpiW^2^) if P > Since j'' x 

b™ and x 6" are both of much larger spatial scale (~ 
size dynamo), they are, from their point of view, either 
evaluated at virtually the same location or the correlation 
function is zero. It follows that 



M 



kmin 



d=^7 



(3'=x5"),(3^x6");a,,(r) , (37) 



where the tensor apg is given by 



'^pqir) = / dV / dT{up{r,t)uq{r + p,t-T)) . (38) 

The integration over p extends over a small radius p ~ 
Ac <C size dynamo, centred on r. Relations ([57)1 and 
psp are a useful starting point for model dynamos with 
small scale turbulence, and we refer to Appendix [C] for 
an application. But often, for example in numerical geo- 
dynamo models, there is no separation of spatial scales. 
The theory developed here remains valid, but the only 
option for computing M*""^" seems to be through mea- 
surements of the turbulent flow u. These permit com- 
putation of time series C^'^{t) from the definition 
after which M'^™^" may be found as an integral over the 
cross correlation function (C'^"'(t)C^"*(t — r)). 



7. MEAN OVERTONE EXCITATION LEVEL 

The idea of determining mean square mode coefficients 
{a^a^*) of a statistically steady dynamo by solving the 
eigenvalue problem corresponding to Eq. p4p may be 
conceptually helpful, but in practice it leads to an ex- 
tremely complicated problem because there are so many 
indices: each single index in Eq. ([34]) actually comprises 
three other indices corresponding to the three spatial de- 
grees of freedom. Approximation methods seem to be 
the only way out. For example, if the geodynamo has 
a dominant dipole mode we may hope that the sums on 
the right hand side of ([M]) are dominated by the term 
m — n — Q. Assuming a statistically steady state we 
then arrive at 



A/c -I- A^ 



,0|2\ 



(39) 



This relation gives the magnitude of the cross correlation 
{a^a^*) relative to the mean square amplitude of the fun- 
damental mode. The mean square overtone amplitude is: 



Since M^^^° 



y u -j^ xb'^ d^r 



/o°°dr(C*'-0(t)C'=o*(t-r)) 
Tc we obtain 



(40) 



-,k\2\ 



^0|2\ 



(41) 



The physical picture behind (HD) is as follows, see Fig.[21 
The turbulent convection u creates new small scale fields 
(overtones) by continuously deforming the dominant fun- 
damental mode. This is encoded in the numerator of 
(I4ip . These overtones are subsequently damped by resis- 
tive effects (5RAfe in the denominator) . The balance of the 
two processes determines the excitation level. Excitation 
of overtone k is effective at locations where f X 6° is 
large, in particular when the turbulent flow u is parallel 
to it. Hoyng and Van Geffen [30] have derived a relation 
similar to (|4ip for a simple dynamo, and these authors 
found that it agrees accurately with numerical results. 

Further progress requires evaluation of (|4ip or of 
j^jkoko^ for which there are basically two options. For 
a model dynamo it may be done analytically. Otherwise 
we need to measure the flow u and compute (|4ip as an 
integral over the correlation function of C^^{t). 

Here we shall restrict ourselves to an estimate of the 
order of magnitude of (j4ip for a spherical dynamo with 
radius R. The volume integration in (j4ip may be esti- 
mated as in (|20p . except that we now do not eliminate 
the number of convective cells N: 



Jv 



(42) 
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so that I / • 



,N-^{k/Rf 



{P/R'^)N ^k^, where (3 is the turbulent diffusion coef- 
ficient. Due to the resistive term in flSl) . JRAfe scales 



approximately as k , so that A; /(— SRAjt) ^ constant 
'-^ B? 1 13. Hence we predict that to order of magnitude 



(|aO| 



/? 1 fc2 

^ TV -SRAfe 



1 



for k <ti 



R 



(43) 



The meaning of k above and in ([^D|) and (PT|) is not im- 
mediately clear as it comprises three quantum numbers n 
(radial) and £, m for the two angular co-ordinates. Since 
k appears when R/k is taken to be the spatial scale of 
6*^, one should assign k = 1 to the fundamental mode. 
For overtones one may think of /c as a geometrical mean, 
k = {n£my/^, but an exact interpretation is of course 
not possible. 

The implication for the geodynamo is that the r.m.s. 
mode amplitude relative to the fundamental mode would 
be of the order of N^^^^ ^0.1, approximately indepen- 
dent of mode number. This does not seem to be an out- 
rageous number, but it cannot be readily compared with 
the data [S^l because these do not distinguish between 
modes of different radial order. 

Estimate assumes effectively that k ^ R/Xc, i-e. 
that R/k, the spatial scale of mode b'^, is much larger 
than the correlation length Ac- For high-order overtones 
with k ^ R/Xc (spatial scale of mode b'^ smaller than 
Ac) we obtain: 



N [kXr 



0|2\ 



R 

for fc > — . (44) 
Ac 



For very high order modes the excitation level approaches 
zero. 



8. DISCUSSION 

We have analysed the statistical properties of the mag- 
netic field generated by a turbulent dynamo in a nonlin- 
early saturated state. The properties of the mean and 
convective flow in this saturated state are supposed to be 
given, which allows us to use linear theory. Starting from 
an expansion in a set of base functions, we have derived 
statistical properties of the expansion coefficients a''{t), 
viz. the means (a*^), the cross correlations (a'^a^*), and 
the autocorrelation functions {a''{t)a^*{t + T)). The con- 
vection may have any spatial scale distribution. Condi- 
tions for validity are (1) a short correlation time, C^^Tc ^ 
1, and (2) the effect of the mean flow v in one correlation 
time is small, S^^Tc <SC 1. The second assumption is made 
for convenience, and the theory can still be deployed if it 
does not hold. These two conditions have been worked 
out for the geodynamo in (|20p and ([2T|) . and seem to 
be well satisfied for the lower multipole coefficients . 
The two main tools enabling our analysis are (1) the fact 
that currents and vector potentials form a biorthogonal 



set in the volume V of the dynamo, and (2) the theory 
of stochastic differential equations. 

Any set of magnetic fields b^{r) may be used for the 
expansion. However, since the theoretical results contain 
averages over the turbulence that also occur in mean field 
theory, it turns out that the eigenfunctions of the dynamo 
equation are a preferred set, in terms of which our results 
assume their simplest form. This is an important point, 
and the reader might easily get the wrong impression, as 
we pay quite some attention to the dynamo coefficients 
aij and /3yfc and to eigenfunctions and eigenvalues of the 
dynamo equation. But this is only done in the interest 
of determining the preferred basis and its properties, not 
because we want to focus our study on the mean field 
concept {B). 



8.1. Nature of the dynamo field 

The physical picture that emerges from our study is 
that the dynamo field is a superposition of transiently 
excited eigenmodes whose coherence time and frequency 
are determined by the corresponding eigenvalue of the 
dynamo equation, as had been surmised by Hoyng [26j. 
For the theory to make sense, all these eigenvalues should 
have negative real parts. This property could be proven 
in the restricted case of locally isotropic turbulence, but 
a general proof is still lacking. The r.m.s. excitation level 
of the modes is given by an equation that determines the 
cross correlations (a'^a^*) up to an overall constant. It 
follows that the relative excitation level of the modes is 
determined by the statistics of the flow and by linear the- 
ory. This is an example showing that linear theory is not 
yet fully exhausted. Determination of the absolute exci- 
tation levels requires explicit inclusion of nonlinear effects 
in the theory. Unfortunately, this equation for {a^a^*) is 
rather complicated and cannot be solved in general. An 
approximate expression for (ja'^p) / (|a°p) exists when 
the fundamental mode is dominant in magnitude, which 
may hopefully apply in the case of numerical geodynamo 
models. 



8.2. Applications 

Application of this work is restricted to kinematically 
stable dynamos, as defined in Sec. [H that are in a quasi- 
steady saturated state. Otherwise, there is no restric- 
tion, and numerical models and laboratory dynamo ex- 
periments are equally well eligible. And although current 
computer technology does not yet permit construction of 
numerical solar dynamo models with adequate resolution, 
such models do also qualify once they become available 
(and turn out to be kinematically stable, which in view 
of their much higher magnetic Reynolds number could 
be a problem) . What is needed is the possibility to mea- 
sure the flow, so that the dynamo coefficients and 
Pijk and thence the preferred basis may be determined. 
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The theory developed here may be extended in several 
directions. One is the computation of the mean rever- 
sal rate of numerical geodynamo models, see Sec. 18.31 
Another would be the statistical distribution of the ex- 
pansion cocfhcients. We did not consider this topic here, 
but it might be of interest as a theoretical underpinning 
of the Giant Gaussian Process approach to geomagnetic 
field modeUing 32, 33]. 



8.3. Mean reversal rate and phase memory 

The geomagnetic dipole reverses its direction at ran- 
dom moments, on average once every few 10^ yr [3l| . 
and numerical geodynamo models exhibit a similar be- 
haviour. In this paper we have argued that the imagi- 
nary part of Afc is the frequency and — l/5RAfc the coher- 
ence time of mode k. We apply that to the fundamental 
mode, and identify the coherence time of the fundamental 
dipole mode with the mean time between reversals. So 
we anticipate that Ao is real (zero frequency) and slightly 
negative (a relatively long coherence time). This leads to 
a simple method for computing the mean reversal rate of 
a geodynamo model: determine the a and P tensors (by 
measuring the flow) and compute Aq. 

Unfortunately, there is a snag. In units of the inverse 
diffusion time P/B? , we will have < — Ao ^ 1, while for 
the overtones — SRAfe > 1. It will be difficult to determine 
the dynamo coefficients accurately enough to attain a 
precision in the eigenvalues of a small fraction of P/B?. 
And this is necessary since the mean reversal rate of the 
geodynamo is ^ 0.01/3/i?^. So this does not seem to be 
a practical way to compute these quantities. However, if 
the fundamental mode is dominant we may use (|39[) with 
/c = ^ = 0, to obtain: 



mean reversal rate — ~ Xn 



M' 



0000 



Here we used to estimate M°°°° - |C™|?.„., Tc - 
{P/B?)/N. For the geodynamo ([45]) is of the right order 
of magnitude, ~ 8 x 10"^'* s~^ or once per 4 x 10^ yr, tak- 
ing for P ~ 100 77 ^ lOOm^s"^ and N ~ 100 convection 
cells. We mention in passing that the determination of 
the mean reversal rate from data is complicated by issues 
such as the time resolution of the data |3J] . 

A related application would be the solar dynamo, 
which has a periodic fundamental mode, 5Ao — 
(22 yr)"^, and according to the theory developed here, 
— l/5RAo would be the coherence time of that mode, i.e. 
the time over which the solar cycle remembers its phase. 
From observations we know that the variability of the pe- 
riod P of the solar cycle ((5P)r.m.s./P ^ 3?Ao/9Ao - 0.1 
[ssl ]. In principle it should now be possible to compute 
this number from theory. 



8.4. The FOSA enigma 

Although mean field dynamo theory is not our main 
focus, we point out that we are now in a position to shed 
new light on the problem of the First Order Smoothing 
Approximation (FOSA). We illustrate the point for the 
solar dynamo, but the argument most likely holds for 
any turbulent dynamo. The FOSA enigma refers to the 
derivation of the dynamo equation, where one is forced 
to make an unjustified approximation, the First Order 
Smoothing Approximation (FOSA). But in spite of this 
the dynamo equation produces very convincing results, 
such as a periodic solar dynamo with migrating dynamo 
waves, a butterfly diagram, etc. In short, mean field 
theory seems to perform much better than one may rea- 
sonably expect, and the question is why? 

A fast road to mean field theory is to interpret pT|) as 
the induction equation, i.e. we put x^B^A — \Jy,v x 
-r^V^ and F{t) = Vxm x , see AppendixIBTDor Ref. 
The result is the dynamo equation ((25|) . and since F = 
(V X Mx)i-.m.s. ~ Wr.m.s./Ac, the short correlation time 
requirement Ftc <€. 1 leads to Wr.m.s.Tc/Ac <C 1 (FOSA), 
which is unlikely to be satisfied in actual dynamos. 

One should be aware of the fact that the condition 
of a short correlation time, Ftc <C 1, is qualitative. It 
is not known by how much Ftc should be smaller than 
unity to avoid that higher-than-second-order correlations 
become important, and this may also differ from applica- 
tion to application. It is possible that in the dynamo case 
Ftc — fur.m.s.Tc/^c with a numerical factor / that is ac- 
tually much less than unity. This could explain why the 
solutions of the dynamo equation behave as if Ftc ^ 1 
even though Ui..„i.s.''c/Ac ^ 1. Unfortunately, it is not 
obvious how this idea may be verified. 

The theory developed here does not consider an aver- 
age of B, but rather averages of expansion coefficients 
of B, and that makes a big difference. This allows us 
to shed some light on the FOSA issue from a different 
perspective. The short correlation time condition (|^ 
translated to the solar dynamo is now 



( A? 



-1/2 



1 



A. 



(46) 



For the number of convection cells we took AirR^d/X'^ 
with d = thickness of the dynamo layer, and %/ Rd — spa- 
tial scale of the fundamental mode, ^/Rd/k for overtone 
k. The remarkable fact is that (|46)) is amply satisfied, 
while FOSA is not (wr.m.sTc ~ Ac). It follows that the 
results of the present paper hold for the solar dynamo. 
The magnetic field consists therefore of a superposition 
of transiently excited eigenmodes of the dynamo equa- 
tion (j25|) with aij and Pijk given by (IB2|1 even though 
FOSA is violated. The FOSA problem would be fully 
solved if we can show that the fundamental mode has a 
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dominant amplitude, that is, in the stationary solution 
{a^a^*) of Eq. ([M]), the largest element should he {\a°\'^). 
We have thus eased the FOSA problem by mapping it 
onto another problem that seems more amenable to a 
quantitative treatment. 

8.5. Outlook 

We are currently in the process of testing the theory 
developed here with the help of a numerical geodynamo 
model. The averages we have computed above depend 
only on a few properties of the turbulent convection: the 
mean flow v and the dynamo coefficients aij and j3ijk- 
These may be inferred from sufficiently long measure- 
ments of the fiow u{t). This allows computation of aij 
and Pijk from their defining relation (jB2[) . and the eigen- 
functions and eigenvalues of the dynamo equation. Once 
we know the basis functions, we may obtain time series of 
C'^Ht) and, finally, M'^"'^" from It seems unlikely, 

incidentally, that M'""^" can be computed with ([37]) and 
because there appears to be no clear separation of 
scales in numerical models. At the same time the dynamo 
field B is to be measured and projected onto [or rather 
the vector potential A is to be projected on j^^ as in (fTTjl ] 
to obtain the expansion coefficients (t) and their statis- 
tical properties. These may then be compared with the 
theoretical predictions. The usual practice in numerical 
dynamo models is that one observes and analyses inter- 
action between local structures. A novel aspect of our 
approach is that the physics of the dynamo may now be 
analysed in terms of interaction of global structures. It is 
hoped that in doing so we gain new perspectives on the 
inner workings of numerical dynamo models. 

Acknowledgements 

I am obliged to Drs. D. Schmitt, M. Schrinncr, R. 
Cameron and Prof. G. Barkema for many useful discus- 
sions. The paper has benefitted considerably from the 
constructive comments of the (unknown) referees. Part 
of this work has been done while I was a guest at the Max- 
Planck-Institut fiir Sonnensystemforschung. I thank the 
MPS for its hospitality and support. 

APPENDIX A: THE FORM OF THE MODE 
EQUATIONS 

Consider first an inner product in V, and construct an 
adjoint set 6' of b' so that these are biorthogonal in V: 
Jy U ■ y d^r = S"-^. The equivalent of ^ is: 

/ 6i-62d^r = / 3i-a2d^r-® 6ixa2-d^cr. (Al) 
Jv Jv JdV 

The surface term is in general nonzero, and is no longer 
cancelled by a second surface term .... Hence currents 



and vector potentials no longer constitute a biorthogonal 
set in V. However, we may still proceed and infer from 
([TU| and (P) that 

a''{t) = [ h^ Bd^r 
Jv 

= [ b'^ -V X {(?; + «)- 77V} X BdV , (A2) 
Jv 

which is pS]) . We may get rid of Vx operating on {• • • } 
by integrating by parts which yields (|12p . however at the 
expense of the following surface term added to the right 
hand side: 

- * b'' X {{v + u) - 77V} X B ■ dV , 
Jav 

which is in general nonzero, i.e. the boundary conditions 
do not force it to be zero. The conclusion is that if we 
use an inner product in V there is no way to infer (jl2p - 
we are stuck with 

Next consider an inner product in V + E, as defined in 
dl]). This leads us to Eq. pH)) . and we then hit the prob- 
lem that it is not clear how to proceed with Jg b^ B d^r. 
The induction equation does not hold in E, and there is 
no obvious choice ioi V = v + u and 77 that one can 
make, so that V x (V — r]V ) x B represents B in E. 
Again we would be stuck, except that we now have the 
biorthogonal current - vector potential mechanism ^ at 
our disposal. We integrate (O by parts and arrive at 
Eq. pip . There is no longer any problem with surface 
terms. Finally we proceed from Eq. (fTT|) to as out- 
lined in the main text. 



APPENDIX B: EXPLICIT FORM OF THE 
DYNAMO OPERATOR D 

There are two ways to arrive at Eqs. (|22ll25p . The 
fast way is to bypass ([19]) and to derive (|25|) and ((24)) 
straightaway. We then insert expansion ^ in (pS)) . and 
use the inner product dH) to isolate Eqs. ([211) -(ESI)- This 
method a disadvantage: it is subject to the FOSA con- 
dition and has therefore a restricted applicability. The 
second method is to start from Eq. P^ and to com- 
pute the correlation function appearing in there, which 
leads us to (|22l) and following equations. This method is 
more laborious and also assumes a short correlation time. 
However, since it uses expansion coefficients of B rather 
than B itself, it turns out to have much wider range of 
validity. This validity issue is not considered here but in 
Sec. [5] and m 



1. Computation via the dynamo equation 

We interpret p7p as the induction equation, and iden- 
tify x = B, A = V x v X -r/V^ and F{t) = V x tt x . 
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Application of (jlSp leads to the dynamo equation (1251) 
with 

Db = I (t> - 7?V) + J dr (m* X V X m*-") | x b 

/•OO 

= (t) — 77V) X b + dr 

( w* X [ (b • V)m*-" - • V)6] ) , (Bl) 

since V-6 = 0, and we assume incompressibility, V-m = 0. 
Time arguments appear momentarily as an upper index: 
M* = u{r,t). Relation (|B1[) is equivalent to with 



and 



(B2) 



where Up = p-th vector component of u{r,t). These are 
the familiar a- and /^-tensors from mean field theory, in 
the FOSA approximation, for incompressible convection, 
and without assuming any flow symmetry. Note that 
0nqs = when n — s. Hence 9 out of the 27 components 
of (3nqs are zero. This is no longer the case when resistive 
effects are taken into account. 

To make the connection with Eq. ([22]) insert expansion 
dSl) in (Ull): 9t(a')b* = {a^)V x Db'^ (summation over 
i and £). Uncurling to dt{a')a^ = {a^)Db^ + V-0, left 
multiplication with j*' and use of the inner product ^ 
produces ([22]) and (|23l) . 



2. Computation from (|19[) 

The disadvantage of the above argument is that the 
requirement of a short correlation time Ft^ ^ 1 leads to 
the restriction u^.m.sJc/ <C 1 (FOSA), while there are, 
unfortunately, several indications that Ur.in.s.''c/Ac 1. 
We can do better if we start from Eq. and compute 
the correlation function. The result is the same, but the 
range of validity is much wider, see Sec. [51 and 18.41 

To evaluate the correlation function in ([TO)) it is useful 
to set up a formal approach, which is really an overkill 
here, but not in the next Section as we deal with resistive 
effects. We consider operators in vector space with a 
left and a right vector gate, for example C = m x , and 
a C b — a uxb. Since (C-b)j — {uxb)i — eiskUghk and 
(C • b)i = Cikbk we have Cik = etskUs- But C has also 
a representation in function space, denoted with upper 
indices. This is C^^{t) = Jy y ■ u x 6*''d^r, defined in 
P^ . While C**'' is an 00 X 00 array of functions of time, 
Cik is a 3 X 3 array of functions of r and t. 

For the computation of the infinite internal summation 
over i in (jl9p . we consider a slightly more general ques- 
tion: given two operators P'"^ and Q''^ in function space 



representation, compute {PQY^ = P'^^Q^^ (summation 
over t) . Here the completeness relation ([5]) figures as the 
essential tool: 



f V d^r / f q, b" dV 

V Jv 



dl d2 Pnrra ^Li fp2 Qpr2 ^2 

— I d V j]^ Pnra ^pms ^ s Qpr 

Jv 

= [ d^r j\ P„„ (V X Q • b'')^ 

JV 

= I d^r y • (P • V X Q) • b'' . (B3) 
Jv 

Notation: f^-^ = n-th vector component of j'(t"i), etc. 
In the second line we insert the completeness relation 
Jp2^mi — ^pmsS{l — 2)Ws2, after which the volume in- 
tegral over 2 = r-z may be done with the help of the 
delta-function. In the third line Vg operates on Qpr and 
6j: as these had originally the argument r2- We see that 
Vx appears at the location of the sum over £. Relation 
(jB3[) tells us that the sum over i in the top line is equal 
to an integral of an ordinary vector expression. It may 
be generalised to several internal summations, e.g. 

{PQTy^ = I d^r y • P • V X Q • V X T • 6\ (B4) 
Jv 

where each V operates on everything to its right. 

The summation over i in (|19p is now almost trivial. 
With the help of (|B3| and C* = m* x we find that 

C'''{t)C'\t-T) ^ [ d^r f ■C{t)-V xC{t-T)-b'^ 
Jv 

= I d^r f -u^ xW X u^-^ X b' . (B5) 
Jv 

Now take Jp°° dr ( ) of (jB5p , combine that with S'^^ from 
(jlSp . and we have proven that the operator between { } 
in IH]) is equal to D''^ = jyd^r ■ D - b^ with D given 
by jlll). 



3. Resistive effects 

Sofar we have ignored the exponential operators 
exp(±AT) in and which is tantamount to ig- 
noring all mean flow and resistive effects on the a and 
/? tensors. Here we investigate the influence of resistiv- 
ity, with the help of the operator technique introduced 
above. Instead of pop our starting point is now 



^ke 



dT{C^'"'{t) [exp(5T): 



C-^{t~T)[eM~ST)Y')\{a') ,{m) 
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with S^^ and C"""(t) given by Uni) and JTll). There are 
three internal summations in (|B6p . over m, n and p, but 
as we shall see the exponential operators contain addi- 
tional summations. As a first step we compute 



[Pexp(5T)Q]*'^- = 

= [P(l + .S'T+i. S^St"^ 



= / dV y • P • exp(V X S r) • V X Q • 6*^ 



(B7) 



In the first and second line P, Q and S are oo x oo arrays 
with numbers as entries, and a bullet • indicates an in- 
ternal summation over a double upper index. To arrive 
at the last line we apply relation (|B4p . As explained in 
the previous section, Vx appears at the location of each 
internal summation; P, Q and S are vector operators, 
e.g. ^ — [v ~ riW)x, and a center dot indicates an in- 
ternal summation over a double lower (vector) index, i.e. 
the usual vector product. Since relation (jB7p is easily 
generalised, we are now in a position to compute: 



dr [ {C{t) exp(S'T) C{t - r)) exp(-S'T) 



I ke 



[ dr [ d^rf ■ (C 
Jo Jv 



3„ //-it gVxSr 



V X C*"^) 



-VxSr , f^e 



d^r "j^ 



/•oo 

/ dr (m* X e^"^' • V X it*-^) 

. (B8) 



X e 



The time argument appears again as an upper index; at 
the second = sign (1) the integration order is reversed, 
(2) we use C = it x , and (3) we stipulate zero mean flow 
so that V X S = —77V X Vx = ?7V^, since V x S operates 
exclusively on vectors having zero divergence. Next, we 
combine dUl) with S^^ from ^ to see that Eq. is 
equivalent to Eqs. and with 

D^{v - 77V) X 

P-OO 

+ dr (m* X e^"^' • V x m*-^) x q-^'^^' . (B9) 
"'0 



The difference with (jBip is in the exponential operators 
exp(±TryV^) that embody the effect of a finite resistivity 
on the a- and /3 tensors. They are evolution operators: 
exp(r77V^)b = B{r,T), where B{r,t) is the solution of 
dfB = ■q\/'^B with initial condition B{r, 0) = 6(r). This 
is usually written in terms of the more familiar Green 
function formalism Q. The operator technique has the 
advantage of being very compact. Explicit expressions 
for the a and /? tensors may be extracted from (IB9p with 
the help of a spatial Fourier transformation of the flow 
u. For this and related matters we refer to I2al. 



APPENDIX C: PROOF OF KAt < 

We close with a proof that the eigenvalues of the dy- 
namo equation of a statistically steady dynamo with lo- 
cally isotropic small-scale turbulence all have negative 
real parts. We begin with Eq. p4|) . take k — i and put 
dt{a^a^*) = on account of statistical steadiness: 



,fe|2 



|2) = _ 1 (^/femfcn ^ Mfc«'='"*)(a'"a"*) . (Cl) 



(summation over m and n, not over k). It is essential that 
we restrict attention to kincmatically stable dynamos. In 
that case the induction equation (on which our results 
are based) will generate the quasi-steady saturated field 
of the dynamo, and {a^a^*} will actually be constant. 

Small-scale turbulence allows us to invoke (|57|) . and 
locally isotropic convection is understood to imply that 
the tensor apq is diagonal: 



Jpq 



(r) = (3X1 Sp, , 



(C2) 



with 



i dV ^ dr {u{r, t)-u{r + p,t-T)) . (C3) 



The expression on the right is written as (3X1 because 
it has the dimension of a turbulent diffusion coefficient 
times a correlation volume. We use /3A|? as a generic 
symbol defined by (|C3[) . It may be a function of position. 
In this approximation we have 



M 



kmkn 



d^r (3X^ if X b™) • if X b")* . (C4) 



V 



Since we now also have that M'^"'^™* — M'^™''", it follows 
that 



SRAfc (|a 



fe|2\ 



d^r pXl {f X b") • [t X b")' 



(a^a"*) 



d^r (3X1 {\f X ^1 



(C5) 



In the last line the summation over m and n has been 
performed with the help of relation It follows that 
3fiAfe < 0, provided (3X^ is positive. This will be the 
case for almost all turbulent fiows u{r, t). We stress that 
the proof allows for a finite resistivity. We have tried to 
extend the proof to more general dynamos, and while we 
have not been unsuccesful yet, we are confident that such 
an extension should be possible. 

The physics behind the above proof is not some kind of 
stability analysis, but rather ordinary phase mixing, cf. 
Sec. 15.31 The field of the dynamo may be represented as 
a sum of oscillators, B — a'b\ The mode amplitudes 
a* — cos{bj^t -I- (^') have a variable magnitude and 
a randomly drifting phase t/s*. The average over an en- 
semble of dynamos {B) — ^^(yl* cos(w*f -|- </3*))b' will be 
zero. It follows that the mean field IB) over an ensemble 
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of statistically steady dynamos is zero, as in ([31]) . This in perturbed or driven oscillators whose mean amplitude is 
turn implies < 0. Another pertinent remark in this zero if the mean energy is constant [Tsl. [26|. 
connection is that there are several examples of randomly 
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